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Beschreibung 

Vorrichtung zur Bestimmung und/oder Uberwachung des Volumen- 
und/ oder Massendurchflusses eines Messmediums 

[001] Die Erfindung bezieht sich auf eine Vorrichtung zur Bestimmung und/oder U 

berwachung des Volumen- und/ oder Massendurchflusses eines Mess-mediums, das 
eine Rohrleitung in einer Stromungsrichtung durchfliesst, mit zumindest zwei Ultra- 
schallwandlern, die Ultraschall-Messsignale in die Rohrleitung aussenden und aus der 
Rohrleitung empfangen, und mit einer Regel-/Auswerteeinheit, die den Volumen- und/ 
oder den Massendurchfluss des Messmediums in der Rohrleitung anhand der Laufzeit- 
differenz der Ultraschall-Messsignale in Stromungsrichtung und entgegen der Stro 
mungs-richtung ermittelt. 

[002] Ultraschall-Durchflussmessgerate werden vielfach in der Prozess- und Automatisie- 
rungstechnik eingesetzt. Sie erlauben es, den Volumen- und/oder Massenstrom eines 
Mediums in einer Rohrleitung beruhrungslos zu bestimmen. 

[003] Bei dem zuvor genannten Laufzeit-Differenz-Prinzip wird die unterschiedliche 
Laufzeit von Ultraschall-Messsignalen in Stromungsrichtung und entgegen der Stro 
mungsrichtung des Mediums ermittelt und ausgewertet. Hierzu werden die Ultraschall- 
Messsignale von den Ultraschallwandlern wechselweise in Stromungsrichtung und 
entgegen der Stromungsrichtung des Mediums ausgesendet und vom jeweils anderen 
Ultraschallwandler empfangen. Aus der Laufzeitdifferenz der Ultraschall-Messsignale 
lasst sich die Fliessge-schwindigkeit und damit bei bekanntem Durchmesser der 
Rohrleitung der Volumendurchfluss bzw. bei bekannter oder gemessener Dichte des 
Mediums der Massendurchfluss bestimmen. 

[004] Hinsichtlich der Typen von Messgeraten wird unterschieden zwischen Ultraschall- 

Durchflussmessaufhehmern, die in die Rohrleitung eingesetzt werden, und Clamp-On 
Durchflussmessgeraten, bei denen die Ultraschall-sensoren von aussen an die 
Rohrleitung mittels eines Spannverschlusses angepresst werden. Clamp-On Durch- 
flussmessgerate sind beispielsweise in der EP 0 686 255 Bl, der US-PS 4,484,478 oder 
der US-PS 4,598,593 beschrieben. 

[005] Bei beiden Typen von Ultraschall-Durchflussmessgeraten werden die Ultraschall- 
Messsignale unter einem vorgegebenen Winkel in die Rohrleitung bzw. in das 
Messrohr, in der/ in dem sich das stromende Medium befindet, eingestrahlt und/oder 
empfangen. Um eine optimale Impedanzanpassung zu erreichen, werden die Ul- 
traschall-Messsignale ttber einen Vorlaufkorper bzw. einen Koppelkeil in die 
Rohrleitung eingekoppelt bzw. aus der Rohrleitung ausgekoppelt. Hauptbestandteil 
eines Ultraschallwandlers ist desweiteren zumindest ein piezoelektrisches Element, 
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welches die Ultraschall-Messsignale erzeugt und/oder empfangt. 

[006] Ublicherweise handelt es sich bei den Ultraschall-Messsignalen, die fiir die 

Volumen- bzw. Massendurchflussmessung verwendet werden, urn breitbandige Pulse. 
Es versteht sich von selbst, dass gerade bei kleinen Nennweiten der Rohrleitung bez. 
des Messrohrs der zeitliche Abstand zwischen dem ausgesendeten und dem 
empfangenen Ultraschall-Messsignal relativ klein ist. Urn hier eine ausreichende Auflo 
sung und damit eine zuverlassige Messung durchfuhren zu konnen, wird das 
Messsignal mit einer Abtastrate abgetastet, die einerseits kleiner ist als die Zeitdauer 
die zwischen dem Aussenden und Empfangen eines Ultraschall-Messsignals liegt und 
die andererseits so klein ist, dass innerhalb der Messimpulsdauer mehrere Abtastwerte 
abgetastet werden. Die Abtast- bzw. Samplerate ist folghch relativ hoch. Die Ab- 
tastwerte bzw. die abgetasteten Amplitudehwerte des Ultraschall-Messsignals werden 
einem A-/D-Wandler zugefuhrt. Eine RegeWAuswerteeinheit, z.B. eine DSP, 
verwendet die Samplewerte bzw. die Abtastwerte dazu, das empfangene Messsignal 
durch eine kontinuierliche Funktion zu interpolieren bzw. moglichst wirklich- 
keitsgetreu zu rekonstruieren. Mathematisch lasst sich dieser Sachverhalt durch die 
stetige Funktion f(t) = f(n [die codierte matematische Formel ist] T) = [die codierte 
matematische Formel ist] darstellen, wobei n = 1, 2, 3, ... - also eine natiirliche Zahl - 
ist und wobei die Koeffizienten [die codierte matematische Formel ist] die zu den 
Zeitpunkten (n [die codierte matematische Formel ist] T) gemessenen Ampli- 
tudenwerte des Ultraschall-Messsignals darstellen. 

[007] Im einfachsten Fall handelt es sich bei der Funktion urn die sukzessive, lineare 

Verbindung von jeweils zwei aufeinanderfolgenden Samplewerten / Abtastwerten. Da 
dieses Verfahren fiir Messungen mit einer gehobenen Messgenauigkeit im Bereich der 
Ultraschall-Durchflussmessung nicht ausreichend ist, ist es bekannt geworden, die 
Lagrange-Interpolation oder die noch komplexere Interpolation nach Levenberg- 
Markart zur Rekonstruktion des empfangenen Messsignals heranzuziehen. 

[008] Die beste und wohl exakteste theoretische Interpolationsmethode fiir aus Ab- 
tastwerten rekonstruierte Signale besteht fraglos in der Verwendung des Shannon- 
Nyquist Theorems, wonach sich eine begrenzte kontinuierliche Funktion, z.B. ein Ul- 
traschall-Messimpuls, durch eine unendliche Summe von gewichteten Spaltfunktionen 
sin (x) / x darstellen lasst. Die korrekte Formel lautet wie folgt, wobei t als Abkurzung 
fiir die Zeit steht und T die Zeitdauer zwischen zwei Abtastwerten reprasentiert: 

[009] [die codierte matematische Formel ist] 

[010] Problematisch bei der Anwendung dieser Formel ist, dass die Rekonstruktion des 
Messsignals nur dann hundertprozentig korrekt ist, wenn die Anzahl der Abtastwerte 
unendlich ist. Um das Shannon-Nyquist Theorem in der Praxis anwenden zu konnen, 
ist es natiirlich erforderlich, die Anzahl der gesampelten Messwerte nach oben zu be- 
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schrSnken. Hierbei ist darauf zu achten, dass die Anzahl der gesampelten Messwerte 
gross genug ist, um eine ausreichende und hinreichende Rekonstruktion des 
empfangenen Messsignals erreichen zu konnen. Problem hierbei ist eine moglichst 
kurze Rechenzeit bzw. eine hinlanglich grosse Kapazitat des hinter der Regel- 
/Auswerteeinheit stehenden Mikroprozessors. Bei den heute bekannten Auswerte- 
verfahren fiir hoch-genaue Durchflussmessungen mittels Ultraschall ist der Ener- 
giebedarf fiir die Lieferung von ('Nahezu'-) Echtzeitmessungen so gross, dass hier nur 
Vierdraht-Messgerate geeignet scheinen. Der Einsatz von sog. Niedrig-Energiegeraten, 
insbesondere von Zweileiter-Ultraschall-Durchflussmess-geraten, war bislang 
aufgrund des hohen Energiebedarfs nicht moglich. Wie bereits angedeutet, ist der hohe 
Energiebedarf vorrangig durch die grosse Rechnerkapazitat des Mikroprozessors bzw. 
der DSP gefordert. Letzdich ist der hohe Energiebedarf eine Folge der aufwendigen 
Auswerteverfahren, welche fur hochdynamische Messungen - insbesondere im 
Bereich der Echtzeit-Messung - mit hoher Messgenauigkeit erforderlich sind. 

[Oil] Der Erfindung liegt die Aufgabe zugrunde, ein Ultraschall-Durchflussmess-gerat 
mit niedrigem Energieverbrauch vorzuschlagen, 

[012] Die Aufgabe wird dadurch gelost, dass die Regel-/Auswerteeinheit mehrere Ab- 
tastwerte ( [die codierte matematische Formel ist] mit i = 1, 2, 3, . . .) eines 
empfangenen Messsignals in definierten Zeitpunkten (t) eines vorgegebenen 
Zeitbereichs ermittelt, 

[013] dass die RegeWAuswerteeinheit den vorgegebenen Zeitbereich des Messsignals 
durch eine stetige Funktion (f(t)) interpoliert, wobei die stetige Funktion (f(t)) durch 
< eine Summe einer vorgegebenen Anzahl (n [die codierte matematische Formel ist] 
N) von Wavelets (W) gebildet ist und wobei jedes Wavelet (W) dem Produkt eines Ab- 
tastwerts mit einer Spaltfunktion ( [die codierte matematische Formel ist] ) mit einer 
Gauss* schen Glockenkurve ( [die codierte matematische Formel ist] , [die codierte 
matematische Formel ist] [die codierte matematische Formel ist] R) entspricht. 
Wie bereits an vorhergehender Stelle erlSutert, besteht das Problem mit der Funktion 
sin (x) / x darin, dass sie fiir praktische Anwendungen viel zu langsam gegen Null 
konvergiert. Die Funktion sin(x)/x ist quasi erst bei Minus-Unendlich bzw. bei Plus- 
Unendlich gleich Null. Durch die Produktbildung mit der schnell abfallenden zuvor- 
genannten Gauss'schen Glockenkurve kann dieses Problem elegant umschifft werden, 
wodurch die Rechenkapazitat der Regel-/Auswerteeinheit ebenso wie deren Ener- 
giebedarf erheblich reduziert werden kann. Die Einsparung liegt im Hinblick auf die 
bislang bekannt gewordenen Verfahren bei einem Faktor, der sich in der Gro 
ssenordnung von 10 bis 100 bewegt. 

[014] Gemass einer vorteilhaften Weiterbildung der erfindungsgemassen Vorrichtung ist 
vorgesehen, dass die RegeWAuswerteeinheit zwischen den Abtastwerten zumindest 
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einen zusatzlichen Abtastwert bestimmt und diesen Abtastwerte bzw. diese Ab- 
tastwerte durch die stetige Funktion approximiert, wobei die stetige Funktion durch die 
Summe einer vorgegebenen Anzahl (n [die codierte matematische Formel ist] N) 
von Wavelets (W) gebildet ist und wobei jedes Wavelet (W) dem Produkt eines Ab- 
tastwerts mit einer Spaltfunktion ( [die codierte matematische Formel ist] ) mit einer 
Gauss'schen Glockenkurve ( [die codierte matematische Formel ist] , [die codierte 
matematische Formel ist] [die codierte matematische Formel ist] R) entspricht. 
Dieses Verfahren ist in anderem Zusammenhang bereits bekannt und tragt den Namen: 
Oversampling. Durch die Zwischenabtastung und Zwischenwertberechnung lasst sich 
eine bessere Auflosung des abgetasteten empfangenen Messsignals und damit eine ho 
here Messgenauigkeit bei der Bestimmung des Volumen- oder Massen-durchflusses 
erreichen. 

[015] Weiterhin schlagt eine giinstige Ausgestaltung der erfindungsgemassen Vorrichtung 
vor, dass die RegeWAuswerteeinheit einen Abszissenwert (t) bestimmt, bei dem ein 
Ordinatenwert der stetigen Funktion (f(t)) einen vorgegebenen Grenzwert erreicht. 
Bevorzugt handelt es sich bei dem vorgegebenen Grenzwert der stetigen Funktion 
(f(t)) urn ein Maximum. Es kann sich jedoch ebenso um einen Nullpunkt, ein Min 
imum oder um einen Wendepunkt handeln. 

[016] 

[017] Bevorzugt ermittelt die RegeWAuswerteeinheit anhand der ersten Ableitung f (t) 
der stetigen Funktion f(t) den Abszissenwert (tmax, tmin) im Maximum und/oder im 
Minimum. Dieser ermittelte Abszissenwert ist nachfolgend jeweils der direkte 
Bezugswert fur die Laufzeit des Messsignals in Stromungsrichtung bzw. entgegen der 
Stromungsrichtung. 

[018] Eine besonders vorteilhafte Ausgestaltung der erfindungsgemassen Vorrichtung 
sieht vor, dass die RegeWAuswerteeinheit den Abszissenwert (tmax), bei dem die 
stetige Funktion ein Maximum erreicht, durch eine lineare Interpolation der ersten 
Ableitung der stetigen Funktion (f(t)) nach folgender Formel erhalt, und wobei tO den 
Abszissenwert einer ersten Schatzung kennzeichnet, bei dem im Zeitintervall (tO - T, tO 
+ T) ein Maximum oder Minimum gemessen wird, und wobei f '(t) die zweite 
Ableitung der stetigen Funktion (f(t)) reprasentiert. Mathematisch lasst sich dieser 
Sachverhalt durch die folgende Formel ausdriicken: 

[019] [die codierte matematische Formel ist] 

[020] Im Zusammenhang mit der erfindungsgemassen Vorrichtung wird es als besonders 
giinstig angesehen, wenn die RegeWAuswerteeinheit zwei Ultrasschall-Messsignale in 
zwei Zeitbereichen miteinander koireliert, die entsprechende diskrete Abtastung von 
Korrelationspunkten durch eine stetige Funktion (f(t)) interpoliert und den Ab- 
szissenwert der stetigen Funktion (f(t)) bestimmt, bei dem der Ordinatenwert das 
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Maximum erreicht, wobei der Abszissenwert ein Mass fiir die Zeitverschiebung 
zwischen den in Stxomungsrichtung und entgegen der Stromungsrichtung gesendeten 
und empfangenen Ultraschall-Messsignale ist. Das Ergebnis dieser Auswertung liefert 
somit direkt die Zeitdifferenz zwischen den beiden in unterschiedliche Richtungen aus- 
gesendeten und empfangenen Ultraschall-Messsignalen. 

[021] Die Messgenauigkeit des UltraschaU-Durchflussmessgerats hangt ent-scheidend 

von der korrekten bzw. optimalen Wahl des Koeffizienten [die codierte matematische 
Formel ist] der Gauss'schen Glockenkurve ab. Urn zu gewahrleisten, dass der Ko- 
effizient [die codierte matematische Formel ist] optimal bestimmt ist, wird er gemass 
einer vorteilhaften Ausgestaltung der erfindungsgemassen Vorrichtung in Abha 
ngigkeit von der Anzahl der Messpunkte (MaxSamplei) bestimmt. Insbesondere 
ermittelt die Rechen-/Regeleinheit jeweils einen optimalen Wert fur den Koeffizienten 
[die codierte matematische Formel ist] in Abhangigkeit von der Anzahl der 
Messpunkte (MaxSamplei) mittels eines mathematischen Simulationsprogramms. 

[022] Bevorzugt ist eine Speichereinheit vorgesehen ist, in der jeweils der optimale Wert f 
iir den Koeffizienten [die codierte matematische Formel ist] in Abhangigkeit von der 
Anzahl der Messpunkte (MaxSample) abgespeichert ist. Beispielsweise sind die Werte 
in einer Tabelle abgelegt. Mit dieser Methode lasst sich wiederum Rechenzeit und 
Energie sparen, da der jeweilige Wert bei Bedarf einfach aus der Tabelle entnonunen 
werden kann. In Abhangigkeit von der jeweils geforderten Messgenauigkeit und/oder 
in Abhangigkeit von der momentan zur Verfiigung stehenden Energie lasst sich so stets 
ein optimales Messergebnis erzielen. 

[023] Aufgrund des geringen Energieverbrauchs ist es moglich, das erfindungs-gemasse 
UltraschaU-Durchflussmessgerat als Zweileiter-Durchflussmessgerat auszubilden. 
Zweileiter-Technologie bedeutet, dass die Energieversorgung des Gerats und die 
Messwertiibermittlung an eine entfemte Leitstelle und ggf . die Konfigurierung und Pa- 
rametrisierung des Gerats von der entfernten Leitstelle her uber lediglich zwei 
Leitungen erfolgt. Da die Verdrahtungskosten iiblicherweise einen relativ hohen Anteil 
der Gesamtkosten ausmachen, lassen sich hier beachtiiche Einsparungen erzielen. 
Aufgrund des geringen Energieverbrauchs ist es dariiber hinaus auch moglich, das 
Durchflussmessgerat mit einer internen Energiequelle auszustatten. Die Kom- 
munikation mit einer entfernten Leitstelle kann dann uber Leitungen oder drahtlos, z.B. 
iiber Funk, erfolgen. 

[024] Die Erfindung wird anhand der nachfolgenden Figuren naher erlautert: Es zeigt: 
[025] Fig. 1: eine schematische Darstellung eines Diagramms, bei dem die abgetasteten 

AmpUtudenwerte des Messsignals gegen die Zeit aufgetragen sind (H Stand der 

Technik), 

[026] Fig. 2: eine schematische Darstellung eines Diagramms, an dem die Arbeitsweise 
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der erfindungsgemassen Vorrichtung visualisiert wird und bei dem die abgetasteten 
Amplitudenwerte des Messsignals gegen die Zeit aufgetragen sind und 

[027] Fig. 3: eine schematische Darstellung einer bevorzugten Ausgestaltung der erfin- 
dungsgemassen Vorrichtung. 

[028] In Fig. 1 sind die Amplitudenwerte [die codierte matematische Formel ist] von 
sechs Abtastwerten, die im Zeitabstand T abgetastet wurden, gegen die Zeit 
aufgetragen. Uber die gestrichelte Linie ist die lineare Verbindung zwischen zwei 
Messpunkten dargestellt Die durchgezogene Kurve f(t) kennzeichnet eine Kurve, die 
nach einem der bekannten Rekonstruktions-Algorithmen, z.B. nach Lagrange oder 
nach Levenberg-Markart gefiinden wurde. Nachteil der bekannt gewordenen Methoden 
ist der hohe Rechenaufwand und damit einhergehend der hohe Energieverbrauch der 
Regel-/Auswerteeinheit bzw. des Mikroprozessors oder der DSP. 

[029] Erfindungsgemass ermittelt die RegeWAuswerteeinheit 1 1 mehrere Abtast-werte ( 
[die codierte matematische Formel ist] mit i = 1, 2, 3, ...) eines empfangenen 
Messsignals in definierten Zeitpunkten (t) eines vorgegebenen Zeitbereichs. An- 
schliessend werden die Abtastwerte in dem vorgegebenen Zeitbereich durch eine 
stetige Funktion (f(t)) interpoliert, wobei die stetige Funktion (f(t)) durch eine Summe 
einer vorgegebenen Anzahl (n [die codierte matematische Formel ist] N) von 
Wavelets (W) gebildet ist und wobei jedes Wavelet (W) dem Produkt eines Ab- 
tastwerts mit einer Spaltfunktion ( [die codierte matematische Formel ist] ) mit einer 
Gauss' schen Glockenkurve ( [die codierte matematische Formel ist] , [die codierte 
matematische Formel ist] [die codierte matematische Formel ist] R) entspricht. 

[030] Im Zeitbereich hat die Formel die folgende Form: 

[031] [die codierte matematische Formel ist] 

[032] Mittels der zuvorgenannten Berechnung lasst sich eine effiziente und schnelle In- 

terpolation bei relativ geringer Rechnerkapazitat erzielen. 

[033] Betrachten wir den Fall des Oversampling und nehmen an, das die Oversampling- 
Rate gleich g ist, wobei gj eine ganze Zahl und grosser oder gleich 2 ist. Folglich ist 
jeder Abtastbereich der Zeitdauer T in § Teil-Abtastbereiche der Zeitdauer T/J 
unterteilt. 

[034] Im folgenden konzentrieren wir uns auf den Zeitbereich -T/2 < t < +T/2. Selbstverst 

andlich konnen auch andere Werte von t analog berechnet werden indem der Ab- 
tastbereich entsprechend verschoben wird. Die Formel fur die Funktion fQT/$) lautet f 
ur i im Bereich: -(0/2+1) [die codierte matematische Formel ist] i [die codierte ma- 
tematische Formel ist] +(0/2), wobei i § 0 sein muss: 

[035] [die codierte matematische Formel ist] 

[036] c(i/g,n) konvergiert mit wachsendem n sehr schnell gegen Null, so dass in der 
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Praxis die Naherung gemacht werden kann, dass c(i/g, n) gleich Null ist 5 sobald der 
absolute Wert von n einen vorgegebenen Grenzwert erreicht Dieser Grenzwert hangt 
letztlich von der geforderten bzw. notwendigen Messgenauigkeit ab und wird 
nachfolgend als MaxSamples gekennzeichnet. Es hat sich gezeigt, dass eine gewii 
nschte Messgenauigkeit sich in der Praxis erreichen lasst, wenn der Wert von 
MaxSamples im Bereich von ca. 3 bis 10 liegt 

[037] Damit ergibt sich die folgende Gleichung: 

[038] [die codierte matematische Formel ist] 

[039] Die Koeffizienten c(i/g,n) werden einmal berechnet und dann in einer Tabelle der 
Grosse (2xMaxSamples + 1) x ( g - 1) abgespeichert. Somit lasst sich nahezu einer 
Verarbeitung der Messwerte in Echtzeit erreichen. 

[040] 

Tabelle 1 



n 


i/@ = -0.25 


i/B = 0.25 


i/g = 0.5 


-5 


XI 


X2 


X3 


0 


Yl 


Y2 


Y3 ! 


r 

5 


Zl 


Z2 


Z3 



[041] Wesendich ist es, anhand der in einem Zeitbereich T abgetasteten Amplitudenwerte 

den Zeitpunkt t zu ermitteln, an dem ein relativer maximaler Amplitudenwert (oder 

max 

auch ein relativer minimaler Amplitudenwert) der Funktion f (t) auftritt. Beispielsweise 
kann dieser Wert zwischen zwei aufeinanderfolgenden Abtastwerten liegen. Ist t 

max 

bekannt, so ist es auch moglich das tatsachliche Maximum f(t ) (oder Minimum) der 

max 

Funktion f(t) zu bestimmen. 

[042] Die Vorgehensweise ist die folgende (siehe auch Fig. 2): Zuerst wird ein Maximum 
(oder Minimum) in dem Zeitbereich bestimmt, indem die Abtastwerte miteinander 
verglichen werden. Nehmen wir an, dass a Q dem maximalen Wert (oder dem minimalen 
Wert) entspricht. Zudem wird angenommen, dass mindestens MaxSamples vor und 
hinter a Q vorhanden sind. Fur praktische Anwendungen liegt der Parameter 
MaxSamples im Bereich von 3 -10. 

[043] Bei dem Abszissenwert t ist die erste Ableitung f (t) der Funktion f(t) also gleich 

max 

Null. Die Formel lautet demnach: f (t ) = 0. Somit ist a das Maximum (oder das 

max 0 
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Minimum) von dem Zeitbereich bzw. Zeitintervall [a], wobei gilt: -T < t < T. 

i 4 ' max 

[044] Der Wert lasst sich durch eine lineare Interpolation der ersten Ableitung gemass 

der folgenden Formel finden: 
[045] [die codierte matematische Formel ist] 

[046] Wie bereits gesagt wird erfindungsgemass folgendes Interpolations- Wavelet 
genutzt: 

[047] [die codierte matematische Formel ist] 

[048] Die erste Ableitung g'(x) und die zweite Ableitung g"(X) der Funktion g(x) lauten: 
[049] [die codierte matematische Formel ist] 

[050] [die codierte matematische Formel ist] 

[05 1] Eine Grenzwertberechnung mit x = 0 fiihrt zu g(0) = 1, g'(0) = 0, 
[052] g"(0) = (-2g-l)/3 

[053] Hieraus ergeben sich fiir f(t), f (t), f'(t) folgende mathematischen Gleichungen: 

[054] [die codierte matematische Formel ist] 

[055] [die codierte matematische Formel ist] 

[056] [die codierte matematische Formel ist] 

[057] Folglich gilt fur t = 0 

[058] [die codierte matematische Formel ist] 

[059] [die codierte matematische Formel ist] 

[060] [die codierte matematische Formel ist] 

[06 1] Wird weiterhin angenommen, dass nur die Abtastwerte innerhalb des Sub- 

Zeitbereichs von a bis a relevant sind, dann wird a =0. wenn ffnll > 

-MaxSampIes Max Samples n u U 

MaxSamples ist. Die letzte Gleichung lautet dann: 
[062] [die codierte matematische Formel ist] 

[063] Wenn der Bereich [a] in den Bereich [a -a Q ] vertikal verschoben wird, so lasst sich 

die zuletzt genannte Formel vereinfachen und man erhalt: 
[064] [die codierte matematische Formel ist] 

[065] Um den optimalen Wert fiir den Koeffizienten a zu finden, sollte t = t gesetzt 

max I max 2 

werden, da sich der Abszissenwert des Maximums bei einer translatorischen Ver- 
schiebung nicht dndern sollte, da er gegenober einer Translation invariant sein mosste. 
Allerdings ist in der strengen Theorie t g t^ (ubrigens sogar dann, wenn 
MaxSamples = _oo ist). Grund ist, dass die Funktion, die den verschobenen Bereich [a. - 
a ] interpoliert keine einfache Translation der Funktion ist, die den Bereich [a] in- 

° i 

terpoliert. Daher ist die Annahme, dass t _« t ist, nur giiltig fiir das nachfolgend 

maxl max 2 " **" 

formulierte Residuum: 
[066] [die codierte matematische Formel ist] 

[067] Der absolute Wert von diesem Residuum ist ein Zeichen dafiir, wie gut das gewa 
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hlte Wavelet fur den praktischen Gebrauch ist. Die Analyse des Residuums mittels 
eines Mathematischen Simulationsprogramms (z.B. Mathcad) ermoglicht es, einen 
optimalen Wert fur gj fiir eine vorgegebene Anzahl von (2xMaxSamples) um a Q herum 
zu finden. Optimale Werte fur liegen in der Grossenordnung von 0.01 bis ca. 0.04. 
Wird ein Wert fiir g eingesetzt, der abseits des optimalen Wertes liegt, so ist die Inter- 
polationsfunktion zwischen zwei Abtastwerten nicht 'glatt', bzw. sie hat signifikante 
Harmonische jenseits der Nyquist Grenze, aber geringe Abweichungen vom optimalen 
Wert (z. B. 5%) haben wenig Einfluss auf das Ergebnis. Letzteres zeigt wie gut die 
vorgeschlagene Vorrichtung zur Durchfiihrung der praktischen Messungen geeignet 
ist. Insbesondere beeinflussen Rundungsfehler das Messergebnis minimal. 

[068] In Praxis werden zwei Koeffizienten-Tabellen [c ] and [d ] (0 < i< MaxSamples), 

aufgebaut. Dann ergibt sich folgender Zusammenhang: 

[069] [die codierte matematische Formel ist] mit 

[070] [die codierte matematische Formel ist] und [die codierte matematische Formel 

ist] 

[071] Die Koeffizientenbereich [c] and [d] weden normiert, so dass gilt: c =1. Der 

1 * i 

optimale Wert fur g, der von der Anzahl der Abtastwerte MaxSamples abhangt, wird in 

eine Tabelle eingelesen. Die Koeffizienten werden berechnet und in einer Tabelle ab- 
gespeichert. Im Betrieb bedeutet das fiir den Mikro-prozessor bzw. die Regel- 
/Auswerteeinheit des Ultraschall-Durchflussmess-gerats, dass nur einfache Rechen- 
operationen wie Additionen, 

[072] 2xMaxSamples Multiplikationen und eine Division durchzufiihren sind. 

[073] Wie bereits an vorhergehender Stelle mehrfach erwahnt, wird zur Messung des 
DurchfLusses mittels eines Ultraschall-Durchflussmessgerats ein Laufzeit- 
differenzverfahren angewendet. Ein UltraschalUmpuls wird in Stromungs-richtung 
(Up) des Messmediums in die Rohrleitung bzw. in das Messrohr eingestrahlt, von 
einem Ultraschallwandler empfangen und nachfolgend werden von einem (high-speed) 
schnellen A-/D-Wandler innerhalb eines vorgegebenen Zeitbereichs mehrere Ab- 
tastwerte [upj gesammelt. Dasselbe Signal wird dann entgegen der Stromungsrichtung 
(Down) in die Rohrleitung bzw. in das Messrohr ausgesendet, ebenfalls von einem Ul- 
traschallwandler empfangen und von dem A-/D-Wandler gesampelt. Innerhalb eines 
vorgegebenen Zeitbereichs werden ebenfalls mehrere Abtastwerte [dn ] gesammelt. 
Die Zeitdifferenz zwischen den beiden Messsignalen ist proportional zur Fliessge- 
schwindigkeit des Messmediums in der Rohrleitung. Die beiden Abtastwerte werden 
miteinander korreliert gemass der Formel: 

[074] [die codierte matematische Formel ist] 

[075] Der maximale Wert dieses Bereichs korrespondiert zu der Zeitdifferenz der beiden 
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zuvorgenannten Ultraschall-Messsignale. Fur eine hochgenaue Berechnung des Volu- 
menstsroms ist dieses Korrelationsverfahren jedoch viel zu ungenau. Daher wird auf 
die Korrelationsfunktion das zuvorgenannte erfindungsgemasse Verfahren 
angewendet. 

[076] Fig. 3 zeigt eine schematische Darstellung der erfindungsgemassen Vorrichtung in 
der Ausgestaltung als Inline-Ultraschall-Durchflussmessgerat 1. Das Ultraschall- 
Durchflussmessgerat 1 ermittelt den Volumendurchfluss bzw. den Massendurchfluss 
des Messmediums 4, das in Stromungsrichtung (S bzw. Up) in der Rohrleitung 2 
fliesst, nach der bekannten Laufzeitdifferenz-Methode. 

[077] Wesentliche Komponenten des Inline UltraschaU-Durchflussmesserats 1 sind die 
beiden Ultraschallwandler 5, 6 und die Regel-/Auswerteeinheit 11. Die beiden Ultra- 
schallsensoren 5, 6 sind mittels einer in der Fig. 1 nicht gesondert dargestellten Befe- 
stigungsvorrichtung in einem Abstand L voneinander an der Rohrleitung 2 angebracht. 
Entsprechende Befestigungsvorrichtungen sind aus dem Stand der Technik hinlanglich 
bekannt und werden auch von der Anmelderin angeboten und vertrieben. Die 
Rohrleitung 2 hat einen vorge-gebenen Innendurchmesser di . 

[078] Ein Ultraschallwandler 5; 6 weist als wesentlichen Bestandteil zumindest ein piezo- 
elektrisches Element 9; 10 auf, das die Ultraschall-Messsignale erzeugt und/oder empf 
angt. Die Ultraschall-Messsignale werden jeweils iiber die Koppelelemente 7, 8 der 
beiden Ultraschallwandler 5; 6 in die vom Medium 4 durchstromte Rohrleitung 2 ein- 
gekoppelt bzw. aus der vom Medium durchstromten Rohrleitung 2 ausgekoppelt. Das 
Koppelelement 7, 8 sorgt fiir eine moglichst gute Impedanzanpassung der Ultraschall- 
Messsignale beim Ubergang von einem Medium in das andere. Mit SP ist ubrigens der 
Schallpfad gekennzeichnet, auf dem sich die Ultraschall-Messsignale in der 
Rohrleitung 2 bzw. in dem Messmedium 4 ausbreiten. Im gezeigten Fall handelt es 
sich urn eine sog. Ein-Traversenanordnung, in der die Ultraschall-wandler 5, 6 
angeordnet sind. Eine Traverse kennzeichnet dabei den Teilbereich des Schallpfades 
SP, auf dem ein Ultraschall-Messsignal die Rohrleitung 2 einmal quert. Die Traversen 
konnen je nach Anordnung der Ultraschallwandler 5, 6 und ggf. unter Einfugen eines 
Reflektorelements in den Schallpfad SP diametral oder chordial in der Rohrleitung 
bzw. in dem Messrohr 2 verlaufen. 
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Anspriiche 

[001] 1. Vonichtung zur Bestiinmung und/oder Uberwachung des Volumen- und/ oder 

Massendurchflusses eines Messmediums (4), das eine Rohrleitung (2) in einer 
Stromungsrichtung (s; Up) durchfliesst, mit zumindest zwei Ultraschall- 
wandlern (5, 6), die Ultraschall-Messsignale in die Rohrleitung (2) aussenden 
und empfangen, und mit einer RegelVAuswerteeinheit (11), die den Volumen- 
und/oder den Massendurchfluss des Messmediums (4) in der Rohrleitung (2) 
anhand der Laufzeitdifferenz der Ultraschall-Messsignale in Stromungs-richtung 
(S; Up) und entgegen der Stromungsrichtung (Down) ermittelt, dadurch ge- 
kennzeichnet, dass die RegeWAuswerteeinheit (11) mehrere Abtastwerte (mit i 
= 1, 2, 3, . . .) eines empfangenen Messsignals in definierten Zeitpunkten (t) eines 
vorgegebenen Zeitbereichs ermittelt, dass die RegeWAuswerteeinheit (11) den 
vorgegebenen Zeitbereich des Messsignals durch eine stetige Funktion (f(t)) in- 
terpoliert, wobei die stetige Funktion (f(t)) durch eine Summe einer vor- 
gegebenen Anzahl (n N) von Wavelets (W) gebildet ist und wobei jedes Wavelet 
(W) dem Produkt eines Abtastwerts mit einer Spaltf unktion ( [die codierte ma- 
tematische Formel ist] ) mit einer Gauss'schen Glockenkurve ( [die codierte 
matematische Formel ist] , [die codierte matematische Formel ist] [die 
codierte matematische Formel ist] R) entspricht 

[002] 2. Vorrichtung nach Anspruch 1, dadurch gekennzeichnet, dass die Regel- 

/Auswerteeinheit (11) zwischen den Abtastwerten zumindest einen zusatzlichen 
Abtastwert bestimmt und diesen Abtastwerte bzw. diese Abtastwerte durch die 
stetige Funktion approximiert, wobei die stetige Funktion durch die Summe einer 
vorgegebenen Anzahl (n [die codierte matematische Formel ist] N) von 
Wavelets (W) gebildet ist und wobei jedes Wavelet (W) dem Produkt eines Ab- 
tastwerts mit einer Spaltfunktion ( [die codierte matematische Formel ist] ) mit 
einer Gauss' schen Glockenkurve ( [die codierte matematische Formel ist] , 
[die codierte matematische Formel ist] [die codierte matematische Formel 
ist] R) entspricht. 

[003] 3. Vorrichtung nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass die 

RegeWAuswerteeinheit (1 1) einen Abszissenwert (t) bestimmt, bei dem ein Or- 
dinatenwert der stetigen Funktion (f(t)) einen vorgegebenen Grenzwert erreicht. 

[004] 4. Vorrichtung nach Anspruch 3 dadurch gekennzeichnet, dass es sich bei dem 

vorgegebenen Grenzwert der stetigen Funktion (f (t)) um einen Nullpunkt, ein 
Maximum, ein Minimum oder um einen Wendepunkt handelt 

[005] 5. Vorrichtung nach Anspruch 1, dadurch gekennzeichnet, dass die Regel- 

/Auswerteeinheit (1 1) anhand der ersten Ableitung f (t) der stetigen Funktion f(t) 
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[006] 



[007] 



[008] 



[009] 



den Abszissenwert (tmax, tmin) im Maximum und/oder im Minimum bestimmt. 

6. Vorrichtung nach Anspruch 1, dadurch gekennzeichnet, dass die Regel- 
/Auswerteeinheit (11) den Abszissenwert (tmax), bei dem die stetige Funktion 
ein Maximum erreicht, durch eine lineare Interpolation der ersten Ableitung der 
stetigen Funktion (f(t)) nach folgender Formel erhalt, und wobei tO den Ab- 
szissenwert einer ersten Schatzung kennzeichnet, bei dem im Zeitintervall (tO - 
T, tO + T) ein Maximum Oder Minimum gemessen wird, und wobei f '(t) die 
zweite Ableitung der stetigen Funktion (f (t)) reprasentiert: [die codierte ma- 
teniatische Formel ist] 

7. Vorrichtung nach einem oder mehreren der vorhergehenden Anspriiche, 
dadurch gekennzeichnet, dass die RegeWAuswerteeinheit (1 1) zwei Ul- 
traschall-Messsignale in zwei Zeitbereichen miteinander korreliert, die ent- 
sprechende diskrete Sammlung von Korrelationspunkten durch eine stetige 
Funktion (f (t)) interpoliert und den Abszissenwert der stetigen Funktion (f(t)) 
bestimmt, bei dem der Ordinatenwert den Maximalwert erreicht, wobei der Ab- 
szissenwert ein Mass fur die Zeitverschiebung zwischen den in Stro 
mungsrichtung (S, Up) und entgegen der Stromungsrichtung (Down) gesendeten 
und empfangenen Ultraschall-Messsignale ist. 

8. Vorrichtung nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass die 
Rechen-/Regeleinheit (11) jeweils einen optimalen Wert fur den Koeffizienten ( 
[die codierte matematische Formel ist] ) in AbhSngigkeit von der Anzahl der 
Messpunkte (MaxSamplei) mittels eines mathematischen Simulationsprogramms 
bestimmt. 

9. Vorrichtung nach Anspruch 8, dadurch gekennzeichnet, eine Speichereinheit 
(12) vorgesehen ist, in der jeweils der optimale Wert fur den Koeffizienten ( [die 
codierte matematische Formel ist] ) in Abhangigkeit von der Anzahl der 
Messpunkte (MaxSample) abgespeichert ist. 
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